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1. INTRODUCTION 

The Switched Reluctance Generator (SRG) offers several advantages over other generators types, 
like high power density, mechanical robustness, no windings and permanent magnets on the rotor, high 
efficiency, performance in a broad range of speeds, low manufacturing costs, and high fault tolerance [1-3]. 
This type of machine is viable for many applications with variable speed demands in harsh environments, as 
in the field of wind power generation, aircraft power systems, battery charging and electrical traction [4, 5]. 

Energies were used to adjust the SRG to low and medium speeds in wind applications, reducing the 
total cost to eliminating the gearbox. Some techniques have been suggested in modern studies to avoid the 
effects of load and speed variation of the voltage generation for SRG systems [6, 7]. Several controllers have 
been introduced for SRG to achieve better dynamic performance. Examples for these controllers are 
Proportional Integral (PI) control which is simple in realization to be employed in SRG control [7, 8]. 

PI controller is widely implemented in the production process as the control strategy. Basically, the 
system response, steady-state error, and the system stability will be improved by PI controller. Moreover, the 
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PI parameters are dependent on the system features. Therefore, the proper or optimum PI parameters are 
necessary to achieve the desired performance. Traditional and intelligent adjustments are methods of tuning 
PI parameters. Ziegler and Nichols suggested the conventional PI adaptation to the formula based on 
observation of the sensitivity, amplitude, and natural frequency of systems [9]. 

For these purposes, increasing the capabilities of PI controllers by adding new features is highly 
desirable. Particle Swarm Optimization (PSO), Artificial Bee Colony (ABC), Differential Evolution (DE), 
Teaching Learning Based Optimization (TLBO), Gravitational Search Algorithm (GSA), Bat Algorithm 
(BA), Pattern Search Algorithm (PSA), Firefly Algorithm (FA), Bio-Geography Based Optimization (BBO), 
Ant Colony Optimization (ACO), Cuckoo Search (CS) Algorithm, Whale Optimization Algorithm (WOA), 
Imperialist Competitive Algorithm (ICA) and Genetic Algorithm (GA) are the methods of tuning based on 
heuristic optimization to improve the efficiency of the mentioned controller types. [10, 11]. 

Grasshopper Optimization Algorithm (GOA) is Saremi's population-based single objective 
stochastic and heuristic optimization technique [12], which emulates Grasshopper's behavior in nature, and 
mathematically models it to solve problems of optimization with contentious variables. Tests were conducted 
with different test functions such as unimodal, multimodal, composite and CEC2005, and real structural 
design problems, show that GOA can efficiently resolve many problems of optimization (also those with 
unknown search areas) [12]. As GOA considers a certain problem of optimization as a black-box and doesn't 
need gradient information from the search area, this enables it a highly appropriate optimization technique in 
different areas for any correctly formulated optimization problem [10]. 

Since the nonlinear nature and/or magnitude of a problem does not affect the GOA and WOA, 
Where early convergence usually shows certain global optimization strategies, the best solution is found with 
faster convergence more efficiently. In this study, taking into account these advantages of the GOA 
algorithm, a GOA-based PI (GOA-PI) controllers proposed for SRG voltage control. It should be noted that 
no such analysis has been suggested before in the literature. Some comparative results between the proposed 
GOA-PI controller and both WOA-PI and PSO-PI controllers will be presented in order to confirm the 
robustness and effectiveness of the proposed method. 

This paper introduces the development procedures of GOA-based PI controllers, WOA and PSO 
optimization techniques. The task of that controller is to generate the turn off angle (Ooff) of the 
magnetization stage of SRG to regulate the generated voltage under different operating conditions, such as 
load and wind speed variations. 


2. OPERATION OF THREE-PHASE SWITCHED RELUCTANCE GENERATOR 

A SRG is a machine that a doubly salient pole supplied by unipolar power converters. The 
configuration of a 3-phase machine with 12 poles on the stator and 8 poles on the rotor is displayed in Figure 
l(a). The asymmetric half-bridge converter (AHB) for a three-phase SRG shown in Figure 1 (b), mainly 
because it enables the machine to be driven both as a generator and as a motor. 
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Figure 1. (a) Machine structure, (b) Asymmetric half bridge converter for a 3-ph SRG 


The windings of the stator are of a concentrated type and simple shape, the rotor has no winding, no 
magnets and low inertia [3, 13]. The characteristics of the SRG depend on numerous features, mainly: 
machine structure (number of phases, number of stator and rotor poles, stator and rotor arcs), magnetization 
characteristic of the laminations, configuration of the converter and methodology of controller [13, 14]. 
Concentrated stator windings are divided into four diametrically symmetrical pairs linked in series to form 


Int J Pow Elec & Dri Syst Vol. 11, No. 1, Mar 2020 : 75 — 85 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 g 77 


phases. Due to The variation in the air gap and non-linear iron magnetization, the machine flux linkage is a 
nonlinear function of the stator current and the angular position of the rotor as in 


p= PCs, 6) (1) 


Where w is the flux linked by the winding, i, is the phase current and @ is the position of the rotor related to 
the aligned position (0 = 0°). 


By the use the finite-element field calculation such as MotorSolve software with defined parameters, 
this complex feature can be realized [15]. The SRM characteristic of magnetization considered in the present 
paper is illustrated in Figure 2(a). The characteristic of magnetization comprises is a family of curves 
signifying the machine flux linkage as a function of current of the phase for various rotor positions from 
aligned position (0 = 0°) to unaligned position (0 = 22.5°). The area under the curve is the co-energy magnetic 
field W., which the relationship could be described, 


W, = f VG, is)di 2) 


The electromagnetic torque T is proportional to the change of the machine's magnetic co-energy We 
at every phase of the SRM, while the machine rotates: 


r = We (3) 


30 ts=constant 





Where L, is the unsaturated phase inductance, then the flux y will be: 


y = L(6).i, (4) 
Then the familiar simplified relationship for SRM torque T is [16] 
_. 1992 aL(@) 
T= a aa (5) 


Where S is positive for motoring and negative for generating modes as shown in Figure 2(b). By 


injecting current into phase windings during the period when the developed torque is negative, electrical 
energy can be generated as 0 changes from 0 to 22.5° as illustrated in Figure 2(b). 
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Figure 2. Magnetization and torque profile of 12/8 SRG 


3. MATHEMATICAL MODEL OF THE SRG 
A MATLAB/SIMULINK tools are used to represent the nonlinear model of 12/8 SRG. The SRG is 
a non-linear electrical system accompanied by a mechanical system which displays the mechanical dynamics 
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of the generator. A single phase terminal voltage V of SRG is linked to the flux linked of the winding as seen 
below [16] 


I dw 
V = lekt; F a (6) 
where, R; is the resistance of the stator. The magnetic flux linkage with the windings is deliberated by. 
t , 
W(t) = fa (V —isRs)dt (7) 
The currents formed by the nonlinear function i,(y,0), are presented as a lookup table ITBL. The 
nonlinear function 7.(i;,8), which is also used as a TTBL look-up table provides for the electromagnetic 


torques created by the stator phases. The two lookup tables ITBL and TTBL used in modeling of the SRG are 
illustrated in Figure 3. 


Te = D5, T Gris) (8) 


Where Teis the electrical torque and j is the number of phases of SRG. The resulting torque equation 
for the SRG can be represented as: 


Te = Tm + J + Bo (9) 


Where Tm is the prime mover torque, J is the inertia of the rotor, and B is the coefficient of friction. 
The summation of electric output power of each phase in one electric cycle Pout is the average power of SRG 
phases 


1 NS ae 
Pot = N Í ijt (10) 
J= 


Where T, Ns, V; and iy are the conduction period of one phase, the number of motor phases, voltage 
and current of Phase /. 
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Figure 3. SRG Lookup tables ITBL and TTBL 


4. PERFORMANCE STUDY OF SRG 

A SIMULINK diagram that used to determine the generator performance is shown in the Figure 4. 
A three-phase asymmetric power converter supplies the SRG. With this structure, the currents of the phases 
can be controlled independently. To implement turn-on and turn-off angles of each phase perfectly, a rotor 
position sensor is used. To calculate the different variables of the system and mark the measurements, 
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measuring devices and signal processing blocks are applied to the system. Table 1 provides the SRG 
parameters used in this paper. 
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Figure 4. Diagram of the simulink for the 12/8 SRG 


Figure 5 shows waveforms of the flux, current, voltage and torque waveforms for turn-on and turn- 
off angles (a =-7°, B = 9.5°). The main simulation results obtained for wind speed 12 (m/s) or N = 1000 rpm. 


Table 1. SRG 
Parameter SRG 
Voltage (V) 600 
Torque (N.m) 650 
Output Power (kW) 100 
Stator poles 12 
Rotor poles 8 
Base speed (r.min”') 1200 
Resistance/phase (Q) 0.03109 


Moment of inertia/(kg.m*) 0.05 
Friction coefficient/(N.m.s) 0.02 


Three-Phase 12/8 SRG ( N= 1000 rpm, Load = 100 KW) 
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Figure 5. SRG flux, current, voltage and torque waveforms at a = -70, B = 9.5oand a constant load 
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The performance of the 12/8 SRG for rotor speed of 1000 rpm with a = -7°, B = 9.5° at a constant 
load of 100 kW is illustrated in Table 2. 


Table 2. Performance of SRG under loading conditions 


Results a= -7°, B=9.5° 
RMS phase current Ims (A) 173.1 
Average phase current Iavg(A) 112.3 

DC bus output Voltage V.(V) 350 

DC bus output Current Iac(A) 286.5 

Output power Pou(kW) 100.2 
Electromagnetic Torque (N.m) 950 


5. CLOSED LOOP OPERATION OF SRG 

A control technique is designed for the closed-loop operation required to control the voltage 
produced. The magnetization period width of the SRG phase varies according to this approach, Maintaining 
the value of 4, fixed (-7°) and controlling the value of the 6,7 through PI controller. If the speed of wind 
turbine or the loads varies, the value of 6, will be changes by the PI controller [ 7, 16, 17 ]. As a control 
strategy in the production process, the PI controller is commonly used. Basically, The response speed, the 
steady-state error and system performance are enhanced by the PI controller. Nevertheless, the PI parameter 
setting is referred to as system process characteristics. The PI controller's transfer function is. 


Gc(S) = Kp += (11) 


The Integral of Time-Multiplied Square Error (ITSE) is used in this paper to evaluate and design the 
proposed controller for the performance index, and it's as given [18]. 


ITSE =f, t [e(t)]? dt (12) 


The PI controller will be tuned off line using both of PSO, WOA and GOA as presented in Figure 6. 
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Figure 6. Structure of SRG voltage controller using PI with evolutionary tuning 


6. MODELLING OF PSO, WOA AND GOA 
6.1. Particle swarm optimization (PSO) 

The PSO-PI controller is recommended in this section. PSO will study the method of tuning PI 
controller parameters. As the following equation, the operating algorithm is based on the best local and 
global solution [18, 19]. 


vitl = wiv + Cyrandx (Pbest; = s) + C rand(gbest = Sk) (13) 

where, vř at k iteration is the velocity of particle i, v‘*tis the particle i updated velocity, w; is the 
weight and different inertia of particle i, Cı and C2 are constants of positive acceleration, S* at iteration k is 
current position of particle i, rand is between 0 and 1 random number, Pbest,; is the i" particle's best previous 


position, and gbest is the best particle in the population among all particles. The new position can, therefore, 
be modified with the present position and updated velocity as [19, 20] 
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The constants positive acceleration C; and C2 are set to 1.6. The weighted inertia w; is set within the 
range (0.4 to 0.9). 


6.2. Whale optimization algorithm (WOA) 

In 2016, Mirjalili and lewis developed an optimization technique that is inspired by the whales’ 
behavior, so-called WOA [21, 22]. Whales are presented as smart and quick animals in finding a prey. A 
whale firstly searches the prey, encircles it and then, by a strategy attacks the prey, called bubble-net hunting. 

WOA assumes that a position is the best solution for a prey and then tries to change their positions 
towards this agent by the other research agents. Like all heuristic-based algorithms, WO simulates the 
behavior of swarm exploration and exploitation. In the beginning of the algorithm, the optimum design 
position is not known in the search space, the initial search agent is expected to be the target prey or close the 
optimum one. The other search agents attempt to change the positions towards the best one when the best 
search agency is determined. The best agent is described below: 


X(t + 1) = X,(t) — A.D (15) 
D = |C.Xp(t) — X(t)| (16) 
Where X p and X stand for best solution of position and the whale current position, respectively, and 


t indicates the number of current iteration. The vectors A and Č are given by equation (17) in terms of a 
random vector 7 in [0, 1] and a shrinking bubbles vector (a) which is decreased linearly from 2 to 0, as 
expressed in (18). 


= 247-4 C=2.7 (17) 
=e (18) 


i ic 





tmax 


The bubble net hunting is also called the spiral system and has a different mechanism, 
because the humpback will travel to the prey in the form of a helix. The spiral motion is expressed 
in Eq (19) 


X(t + 1) = D’.e”.cos(2nl) + Xp(t) (19) 


Where D’ = |X p(t) — X (t)| and the distance between the prey and the ith whale is determined, b is 
the logarithmic spiral constant, and / in [-1, 1] is a random number. The developers of WO assumed that 
probability (p) determines the movement type, so any value of p < 0.5 means that shrinking mechanism 
applies. Then, the spiral mechanism applies in case of p > 0.5. This allows WO to conduct a global search. 
Additionally. Readers may refer to WO for more information. [21]. 


6.3. Grasshopper optimization algorithm (GOA) 

GOA is a newly proposed single target, a population-based heuristic algorithm that emulates 
grasshopper swarms ' behavior in nature and models them mathematically to optimization problems with 
litigating variables [12]. Among Grasshoppers, the algorithm simulates repulsion and attraction forces. While 
forces of repulsion allow grasshopper to search space, forces of attraction urge them to exploit promising 
regions. So that the exploration is balanced and extraction process, GOA was provided with a coefficient 
decreasing the grasshoppers comfort zone. This allows GOA not to get trapped in local ultimate and find a 
reliable global optimum calculation. Because the best solution the swarm has achieved so far is regarded as a 
chased objective, the grasshoppers have a strong opportunity of finding the global optimum through the 
improvement of the target over the course of the iterations [23, 24]. The GOA equation of position update is 
given by 





bq — lb 
Ke = r (jar ese 5( Xf — 


j+i 


ra) FT, 20) 


Where, X a is the position in d-th dimension the current solution, 7 is a coefficient of diminishing, 
which narrows the area of comfort, repulsion and attraction zone, ubg in d-th dimension upper bound, lbg in 
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d-th dimension lower bound, S describes the social forces among grasshoppers, d;jshall be absolute value for 
the distance among j-th grasshopper X; and i-th grasshopper X; and T4 in d-th value of dimension target, 
which solution has been found to date. Eq (20) shows that a grasshopper's next position depends on his 
position current, all other grasshopper’s position, and the target position. in Eq (20) the function of social 
forces is defined as 


s= jer —e4 (21) 


While, f indicates the attraction strength and / is length scale of attractive. In proportion to the 
iterations number to balance of exploration and exploitation, parameter r in Eq (20) should be reduced. It 
facilitates the use as the increasing number of iterations. It also reduces the comfort zone according to the 
iterations number and is measured accordingly [25] 


r= Vmax 


-t Umax — min (22) 
T 

Wherever fmax 18 the value of maximum, fmin 1s the value of minimum, ¢ is the current iteration, and 

T is the number of iterations maximum. In this paper fmax and Tmin 18 1 and 0.00001 respectively. 


7. SIMULATION RESULTS 
7.1. System response with PI controller tuned by PSO, WOL and GOA 

To maintain the output voltage of the SRG at 350V, for the output checking of the controller 
parameters PSO-PI according to the evaluation, the parameters of PSO are used as follows: 
— Size of population = 25; 
— Wax = 0.9, Wmin= 0.1; 
= C {= C t= 1.6; 
— Iteration = 60; 

For the WOA, there are only two control parameters; @ which was directly reduced from 2 to 0 and 
the random vector 7 in [0, 1]. The best results were obtained using the typical values, Iteration =60 

To verify the Performance of the GOA-PI controller parameters, the following parameters of GOA 
are used: 
— Number of search agents = 25 
— Number of iterations =60 

Figure 7. shows the convergence curves of the WOA, PSO, and GOA algorithms. It can be clearly 
seen that GOA converges to the global value faster than the other algorithms for the problem under 
consideration. By heuristic tuning methods for the tuned and optimized PI controller, the parameters obtained 
after the simulation process are illustrated in Table 3. 
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Figure 7. The fitness function versus number of iterationa 


Table 3. Parameters of optimized PI 


Controller Type Kp Ki 
PSO_PI 0.0231 2.0121 
WOA PI 0.0985 2.9537 
GOA PI 0.0665 2.745 
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As shown in Figure 8(a), the proposed controller can hold the voltage generated at the rated value of 
350 V. Figure 8(b) shows the power supplied to the load. The performance features that define the transient 
response of a unit step input are maximum overshoot, settling time, rise time and steady state error. Thus, the 
implications of these features are shown in Table 4. The performance of the proposed GOA-PI controller in 
comparison with other controllers is verified for these effects. 
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Figure 8. Step response of the system with PI controller tuned by PSO, WOA and GOA 


Table 4. Results of the system for transient response analysis. 
Controller Type Overshot % Settling Time (sec) Rising Time (sec) Peak Time (sec) ITSE value 


PSO PI 0.258 0.0569 0.0249 0.0919 0.2107 
WOA PI 0.2388 0.0508 0.022 0.0895 0.20686 
GOA PI 0.238 0.0506 0.0214 0.0869 0.2068 


GOA-PI has better results than PSO-PI for maximum overshoot by 7.75 % and 0.335% compared to 
WOA-PI. For settling time, GOA-PI has better results by 11.07% compared to PSO-PI and 0.394% compared 
to WOA-PI. For rise time, GOA-PI has better results by 14.056% compared to PSO-PI, 2.73% compared to 
WOA-PI. For peak time, GOA-PI has better results by 5.44% compared to PSO-PI and 1.822% compared to 
WOA-PI. In the same Table 4 are also indicated the performance index values for different controllers. The 
proposed controller gives minimum of ITSE value in comparing with the other controllers, as seen in this 
table. This results confirm that the GOA tuned controller has better performance than the other PSO and 
WOA tuned controllers. 


7.2. Testing of the proposed controller 

To clarify the robustness of the suggested controller, two scenarios will be occurred in the load at 
0.3 second. In Figure 9, the load is decreased suddenly from 100 kW to 58 kW. Furthermore, the load is 
increased suddenly from 64 kW to 100 kW, as shown in Figure 10. 
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The system voltage response with optimum PI controller tuned by GOA, WOA and PSO techniques 
using /TSE fitness function are shown in Figures 11 and 12. Due to the load variation, the SRG terminal 
voltage is enhanced closely to 350 V with small overshoot and recovery time. The results obtained verify the 
effectiveness of the GOA-based voltage control system. 
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In addition, the changing of the wind speed is decreased from 12 m/s to 10 m/s at 0.4 second as 
illustrated in Figure 13. During this event, the terminal voltage of the SRG terminal voltage is regulated 
closely to 350 V with small overshoot and recovery time. Also, it is revealed that the system voltage response 
with optimum PI controller tuned by GOA, WOA and PSO techniques using /7SE fitness function, as shown 
in Figure 14. The results obtained verify the effectiveness of the proposed GOA based controller for voltage 
control system. 
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8. CONCLUSION 

This paper presents output voltage control of SRG based on wind turbine of 100 KW with electric 
grid using a PI controller tuned by both GOA, WOA and PSO. For the parameter tuning technique, the 
control principle is to direct the switches of the electrical converter, to provide optimal parameters for the 
ITSE-based PI controller, the GOA algorithm is frequently used. The controller is tested in three cases, 
increases the load, decreases load and changes in wind speed. The simulation results show that the proposed 
GOA-PI controller can find optimum controller parameters quickly and efficiently. The GOA-PI controller is 
also compared with the WOA-PI and PSO-PI controllers using the results of transient analyses and validity 
analysis. The results of simulation using a GOA-PI controller for SRG based wind turbine show better 
performance of voltage control than the other tuning methods. 
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